arXiv:astro-ph/9502092vl 22 Feb 1995 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 


Printed 1 February 2008 


Einstein Observatory evidence for the widespread baryon 
overdensity in clusters of galaxies 

D.A. White and A.C. Fabian 

Institute of Astronomy, Madingley Road, Cambridge CBS OHA 


Accepted 1994 September 29. Received 1994 July 21; in original form 1994 March 28 


ABSTRACT 

We analyse the X-ray surface brightness profiles of 19 moderately distant and lu¬ 
minous clusters of galaxies observed with the Einstein Observatory . Our aim is to 
determine cluster gas masses out to radii between 1 and 3 Mpc, and to confirm the ap¬ 
parent conflict, if Oq = 1, between the current calculations of the mean baryon fraction 
of the Universe expected from standard primordial nucleosynthesis, and the fraction of 
the mass in clusters which is in gas. Our analysis shows that baryon overdensities in 
clusters are much more widespread than only the Coma cluster with which S. White 
& Frenk originally highlighted this problem. The uncertainties involved in our analysis 
and some cosmological implications from our results are briefly discussed. For a refined 
sample of 13 clusters we find that the baryon fraction for the gas within IMpc lies 
between 10 and 22 per cent. 

Key words: galaxies: fundamental parameters, intergalactic medium, observations - 
cosmology: dark matter - X-rays: galaxies. 


1 INTRODUCTION 

Current constraints, from a comparison of standard, homo¬ 
geneous, primordial nucleosynthesis calculations with light- 
element abundances by Walker et al. (1991, see also Olive 
et al. 1990 and Peebles et al. 1991), place tight limits 
on the baryon density parameter of the Universe, fib = 
0.05 ± 0.01/i^g^ {Ho = 50/i5okms“^ Mpc“^). This implies 
a mean baryon density of about < 6 per cent of the critical 
closure density if /150 = 1, which is much lower than the pre¬ 
vious upper limit of fib < O.lQh^g determined by Yang et 
al. (1984). As first pointed out by S. White & Frenk (1991), 
the recent calculation of the mean baryon fraction conflicts 
with the X-ray determinations of the gas fraction of mass in 
clusters if flo = 1. This fraction should equal flb/flo if dark 
matter is distributed similarly to the X-ray emitting gas. 
They, and more recently S. White et al. (1993), noted that 
hot gas contributes approximately 20 per cent to the total 
mass of the Coma cluster within approximately 3 Mpc, indi¬ 
cating that flb/flo ~ 0.3. Baryon fraction estimates for the 
Coma cluster have been extended to 5 Mpc, utilizing the ‘un¬ 
limited’ field of view of the ROSAT All-Sky Survey (Briel, 
Henry, & Bohringer 1992), where the gas mass content is 
then 30 per cent. Thus at face value, either fig ~ 0.2, dark 
matter and baryons have different distributions on cluster 


scales, or the abundance measurements and/or calculations 
of fib are incorrect. This last possibility is unlikely since it 
involves the best-understood physics. 

Over a decade ago, Ku et al. (1983) determined that the 
baryon fraction within 1.9 Mpc of CA0340 — 538 was greater 
than 10 per cent, while Stewart et al. (1984) found varia¬ 
tions between 3 and 20 per cent within the central 0.5 Mpc 
of 36 clusters, and that a significant number of clusters have 
baryon fractions of at least 10 per cent. Stewart et al. also 
noted that the baryon fraction increase with radius, a result 
supported by the analysis of Einstein Observatory data by 
Forman & Jones (1984) which showed that the scale-height 
of the gas distribution in clusters is generally larger than 
that of the gravitational mass. Edge & Stewart (1991b) also 
determined baryon fractions of up to 20 per cent from EX- 
OS AT observations of 36 clusters of galaxies. However, none 
of these studies noted any conflict between the X-ray deter¬ 
minations of the baryon fraction and the constraints from 
standard primordial nucleosynthesis; the calculated limit of 
the baryon fraction at that time was fib < 0.19 /i^q^. 

A recent study of the Shapley supercluster by Fabian 
(1991), where gas mass and luminosity relations from For¬ 
man & Jones (1984) are extrapolated, indicates that the 
baryon fraction there is greater than 18 per cent over a re- 
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gion of 37 Mpc in radius. This over-density of a factor of 3 
implies that the baryons must have been accumulated from a 
region that is at least 40 per cent larger in radius, if flo = 1- 
This then implies that the Shapley region must be bound, or 
has at least retarded the Hubble flow over this region, and 
creates problems for the theory of the formation of large- 
scale structure in currently favoured models. 

If baryon over-densities are common in clusters, as we 
aim to show in this paper, then perhaps the most obvious so¬ 
lution is that flo < 1 (e-5- for baryon fractions of 30 per cent 
fIb/fIo~0.06/0.2 = 0.3 — see also S. White & Frenk 1991, 
and S.White et al. 1993). However, this solution disagrees 
with the strong evidence for flo = 1 from cluster evolution 
and substructure studies (Richstone, Loeb & Turner 1992), 
and the estimates obtained from the POTENT analysis of 
IRAS galaxy density fields and peculiar velocities (Nusser & 
Dekel 1993, Dekel et al. 1993, Dekel & Rees 1994). If flo = 1, 
there are several possible solutions or implications from the 
baryon over-density problem. 

(i) The X-ray emitting gas has been concentrated with 
respect to the dark matter and the total cluster masses are 
much higher, or there is clustering of dark matter on larger 
scales, e.g. in a mixed dark matter Universe. 

(ii) The X-ray determined gas masses are over-estimated, 
e.g. due to clumping of the X-ray gas. 

(iii) The current examples of large baryon over-densities 
in clusters are unrepresentative of clusters in general. 

(iv) The cosmological constant. A, is non-zero and con¬ 
tributes to the density parameter such that flo = ^matter + 
Ua = 1 (e.g. see the review by Carrol, Press & Turner 1992). 

(v) fib can be higher if the Universe is inhomogeneous at 
the time of nucleosynthesis. 

(vi) The new calculations of standard primordial nucle¬ 
osynthesis or the primordial abundance measurements are 
incorrect, or they are very less tightly constrained, allowing 
fib ~ O.Sh^Q (e.g. if some abundance determinations are 
incorrect). 


Recent 

work on inhomogeneous nucleosynthesis (Jedamzik, Fuller, 


& Mathews 1994) now shows that (v) is not viable, and Kvi)| 
seems unlikely given the physics involved is well understood. 
Solution (ii) is unlikely to cause a significant discrepancy, as 


shown later in this paper (see also McHardy et al. 1990). 
We shall discuss solution |^, which implies that there are 
significant masses of gravitating matter outside the regions 
of the X-ray emitting gas. The main focus of this paper is 
to investigate point (iii) by determining gas masses to large 
radii in a number of clusters. Our determinations have been 
made from X-ray image-deprojection analysis of 19 clusters 
of galaxies which were observed with the Einstein Observa¬ 
tory Imaging Proportional Counter (IPC). These are X-ray 
luminous (I/xboi~5 x 10"^"^ergs“^) and moderately distant 
{z > 0.05) clusters which are easily covered by the IPC field 
of view, and have no strong contaminating sources to their 


surface brightness profiles. This enables their gas masses to 
be well determined out to between 1 to 2.5 Mpc. 


2 DEPROJECTION ANALYSIS AND RESULTS 

We have used the X-ray image deprojection technique (pio¬ 
neered by Fabian et al. 1981). This method assumes a spher¬ 
ical geometry for the cluster, and enables the volume count 
emissivity from the hot intracluster gas to be determined 
as a function of radius. The properties of the intracluster 
medium (ICM) can then be determined after corrections for 
attenuation of the cluster emission due to absorption from 
intervening material, and assumptions on the form of the 
underlying gravitational potential of the cluster. Detailed 
descriptions and recent examples of the current analysis 
method and procedure can be found in the description of the 
analysis of RGAAT Position Sensitive Proportional Counter 
(PSPC) and High Resolution Imager (HRI) data, on A478, 
by Allen et al. (1993) and D. White et al. (1994), respec¬ 
tively. 

This analysis differs from previous deprojection analy¬ 
ses, which have concentrated on investigation of the cool¬ 
ing flow properties of clusters, as such studies required rel¬ 
atively small radial bins to resolve the cooling time of the 
intracluster gas at the very centre of the cluster. We are in¬ 
terested in accurate determinations of the total gas masses 
in clusters to large radii and so large radial bins, which im¬ 
prove the signal-to-noise of the data, enable deprojection 
to greater radii. The clusters in our sample were selected 
to be bright and moderately distant so that the emission 
can be followed to sufficient radius within the field of view 
of the IPC. The X-ray emission of the cluster should also 
be relatively symmetrical with no significant contamination 
from sources which will produce significant errors in the gas 
mass determinations. Most of the clusters which we selected 
do appear fairly spherically symmetric and smooth in the 
central regions, although there are some clusters where the 
emission in the outer regions is less regular (notably A1763, 
A3186, A3266 and A3888). However, as these results are not 
significantly different from the clusters which appear very 
regular (such as A85, A478, A644, A1795 and A2009), we 
believe the results provide a good statistical indication of the 
baryon fraction in clusters, despite morphological details. 

With the above criterion we formed our sample of 
19 clusters, and obtained surface-brightness profiles from 
C. Jones, W. Forman and C. Stern at the Harvard- 
Smithsonian Center for Astrophysics. We selected IPC 
rather that HRI data for this analysis due its larger field 
of view and superior quantum efficiency. The poorer spa¬ 
tial resolution of the IPC, as we have noted is inconsequen¬ 
tial. The point-spread function of the IPC, which is approxi¬ 
mately 1 arcmin (Gaussian width), corresponds to relatively 
large radii at the moderate redshifts of the clusters in our 
sample. The data were corrected for the effects of the tele- 
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scope vignetting, and the background contributions were es¬ 
timated from the region just outside the maximum radius 
of each deprojection. This ensures that our estimates of the 
cluster gas masses are conservative as there may be some 
cluster emission outside the selected maximum radius. The 
deprojection of the surface brightness profile of each cluster 
requires the additional information, shown in Table such 
as the column density along the line-of-sight, the X-ray tem¬ 
perature and velocity dispersion of each cluster. 

The total attenuation of X-rays from the cluster is de¬ 
pendent on the absorption (note we use photoelectric ab¬ 
sorption cross-sections given by Morrison & McCammon 
1983) within our Galaxy and intrinsic absorption. D. White 
et al. (1991b) have shown that many clusters appear to have 
intrinsic absorption, but as we do not have such information 
for our whole sample we use estimates for the Galactic con¬ 
tribution taken from the 21cm determinations by Stark et 
al. (1992). The effect of possible excess absorption will be 
shown in Section but we note here that our prescription 
leads to conservative gas mass estimates. 

The Einstein Observatory IPG data do not have suffi¬ 
cient combined spatial and spectral resolution to enable an 
accurate empirical determination of the temperature profile, 
from which the gravitational potentials of the clusters may 
be determined. This means that the deprojection technique, 
which could otherwise be used to directly determine the to¬ 
tal gravitational mass of the cluster as a function of radius, 
actually requires the form of the gravitational potential to be 
specified. The deprojection results are then calibrated using 
the only widely X-ray observed property of the intracluster 
medium — i. e. spatially averaged cluster temperatures from 
broad-beam detectors (Edge & Stewart 1991a and David et 
al. 1993). 


The form of the gravitational potential that we have 
chosen is that of a true isothermal sphere. This produces 
comparatively conservative mass estimates (compared to a 
King-law distribution), and can be parameterised using ob¬ 
servational data, such as the optical velocity dispersion. In 
our standard deprojection model we used a two-component 
true-isothermal potential, each parameterised by a velocity 
dispersion and core-radius, with one potential for the clus¬ 
ter and another for a central cluster galaxy. Only the cluster 
potential was varied; the galaxy potential was fixed with a 
galaxy velocity dispersion of 350 km s“^ and a core-radius of 
2 kpc. The effect of uncertainties in the cluster velocity dis¬ 
persion, the effect of the mass from a central galaxy, and the 
use of different gravitational mass distributions on the re¬ 


sults were all investigated, and are discussed in Section 3.2 


First we discuss the choice of cluster velocity dispersions and 
core radii. 


value from the following relationship between the velocity 
dispersion and observed X-ray temperature: 

a = 376 [Tx (keV)]° ®^® hms”! . (1) 

This relationship was determined (D. White et al. in prepa¬ 
ration) using the ‘orthogonal distance regression’ algorithm 
(see the ODRPACK V2.01 software by Boggs et al. 1990, dis¬ 
cussed in relation to astronomical data analysis by Feigelson 
& Babu 1992), and accounts for errors in both axes of the 
data — an essential feature when the errors in both dimen¬ 
sions are significant. The final velocity dispersions that were 
used, and the source of these values, are given in Table |^. 

Suitable values for cluster core radii are more difficult 
to obtain than velocity dispersion values. Although they are 
available from the X-ray surface brightness profiles of clus¬ 
ters, and are less prone than optical values to contamination 
from sub-structure within the cluster, they can be affected 
by the presence of a cooling flow (which enhances the X- 
ray emission within the central 200 to 300 kpc, as shown 
by Forman & Jones 1984). Therefore we have not used the 
values for core-radii given in the literature, but used the 
core radius as a free parameter because it can significantly 
alter the shape of the gravitational mass distribution. As 
the best and most widely available cluster temperatures for 
most of our sample are only spatially-averaged values for 
the whole cluster, determined from broad-beam detectors, 
we vary the core-radius and outer pressure to produce a 
temperature profile that is as consistent with the observed 
value over as large a radius of the cluster as possible, i.e. 
a flat temperature profile. This tends to overestimate the 
temperature at the centre of a cluster in a cooling flow clus¬ 
ter, but will lead to conservative estimates of the gas mass, 
as Mgas oc Tx ^ 'pjjg gjjaj selections of core radius and 
outer pressure used in each cluster are given in Table ^ 

Note, we do not assign a particular significance to the 
core-radii that we have used in this analysis; it was essen¬ 
tially used as a parameter to obtain flat deprojected temper¬ 
ature profile for each cluster. This also produces conservative 
gas mass estimates, because the temperature at the centre 
will be hotter than expected in a cooling flow cluster. This 
may not represent the true form of the temperature profile, 
and our resulting core-radii may be somewhat misleading. 
This may explain some of the large core-radii, although it 
may also be due to unresolved physical substructure in the 
X-ray emission. We also note that as the baryon fraction 
varies with radius according to the core radius used, as can 
be seen in Fig. We have therefore quoted our results at the 
maximum radius of each deprojection to ensure the results 
are not affected by the core-radii that were used. 


The cluster velocity dispersions were chosen from the 
literature where available. However, when we could find no 
suitable value, or there appeared to be some problem ob¬ 
taining a satisfactory deprojection results, we obtained a 


The baryon fractions that we determine from the depro¬ 
jected value of Mgas /Mgrav at the maximum radii are given 
in Table ^ They do not include the stellar contribution to 
the baryon content (perhaps an extra 5 per cent). The re¬ 
sults indicate that there is a wide variation between approxi- 
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Table 1. Input Data. 


No. 

Cluster 

z 

Galactic 21 cm 

Temperature (keV) 

Gravitational Potential 




Ah (10^1 cm-2 ) 

Reference 

Deprojected 

rcore (Mpc) 

a (km s ^ ) 

1. 

A85 

0.0521§ 

0.30 

^+0.4(0.2) 
'^'^-0.5(0.3) 

(6-in) 

0.10 

749§ 

2. 

A401 

0.0748°° 

1.11 

* +1.1(0.6) 
'•^-0.9(0.6) 

(S'Sla?) 

0.60 

1112'^^ 

3. 

A478 

0.0881§ 

1.36 

tg8+i.i(o.6) 
^•^-1.0(0.6) 

(7.2lJ:«) 

0.20 

904§ 

4. 

A545 

0.1530°° 

1.14 

^•^-2.8(2.1) 

(e-ns) 

0.65 

925'^-^ 

5. 

A644 

0.0704°° 

0.73 

^+3.0(1.1) 
'•^-1.2(0.8) 

(7.24^:^) 

0.40 


6. 

A665 

0.1816°° 

0.42 

^+1.0(0.6) 

®'^-0.8(0.4) 

(9n:8) 

1.00 

1201°° 

7. 

A1413 

0.1427°° 

0.20 

* +0.5(0.3) 

®'^-0.5(0.3) 

(lo.ojjij) 

0.50 

1193'^^ 

8. 

A1650 

0.0840°° 

0.15 

.+2.7(1.3) 
^•^-1.5(1.0) 

(e-itas) 

0.35 

925'^^ 

9. 

A1689 

0.1810°° 

0.19 

no 1+2-7(54) 
^'^•-^-1.5(1.0) 

(10.9t°;2) 

0.40 


10. 

A1763 

0.1870^ 

0.09 

*6 9+'^ 

^'^-3.6(1.9) 

(7.il°;®) 

0.70 

1043'^y^ 

11. 

A1795 

0.0621§ 

0.12 

1 +0.4(0.2) 

^•-^-0.5(0.3) 


0.20 

773§ 

12. 

A2009 

0.1530°° 

0.33 

*7 q + '^(4-4) 
^•^-2.9(2.1) 

(8.0«;5) 

0.40 

1112'^y^ 

13. 

A2029 

0.0765§ 

0.24 

*^.+1.4(0.8) 
^•^-1.0(0.7) 


0.30 


14. 

A2142 

0.0899°° 

0.39 

1 +2.0(1.2) 
-^-*■•^-0.7(0.4) 

(10.411°) 

0.40 

1295* 

15. 

A2163 

0.2030* 

1.10 

q+1-1^-7) 
^■^•^-0.8(0.5) 

(14-2^3) 

0.60 

1509'^>^ 

16. 

A2319 

0.0559* 

0.86 

* +1.4(0.8) 

^'^-1.1(0.7) 

(ll+lo) 

0.60 

126 

17. 

A3186 

0.1270* 

0.60 

*5.9 

(6.7tl«) 

0.50 

960'^y^ 

18. 

A3266 

0.0545'^ 

0.30 

9+0.6(0.5) 
'^•^-0.6(0.4) 

(6.8i;;«) 

0.80 

985'^y^ 

19. 

A3888 

0.1680^^ 

0.11 

*7.9 

(7.+?:o) 

0.50 

1120'^y^ 


This table contains the input data required for the cluster deprojections. The first temperatures given are reference values (with 5th 
and 95 percentile confidence limits and Icr standard deviations in the brackets) obtained from the literature. In the next column are 
the spatially-averaged emission-weighted (0.4 — 4keV) temperatures from the deprojected temperature profiles (these are median values 
with 10th and 90th percentile limits given in brackets). A comparison of these two columns shows the accuracy of the deprojection 
calibration with respect to the reference temperatures. The velocity dispersion values written in italics with the superscript Tx refer to 
values interpolated from the X-ray temperature values using equation stated in the main text. Note, we have used X-ray temperature 
interpolated velocity dispersions for A401, A2009 and A2029, as we were unable to obtain a fiat temperature profile from the literature 
values. The velocity dispersion for A401 was reduced from 1290°° km s“^ , for while A2009 and A2029 the velocity dispersion was 
increased from 804°° kms“^ and 786^ kms“^ . The superscripts refer to: * David et al. (1993); f Edge & Stewart (1991a); | Forman & 
Jones (private communication); § Zabludoff, Huchra & Geller (1990); oo Struble &; Rood (1991); • Quintana & Lawrie (1982); A Arnaud 
et al. (1992); 0 Noonan (1981); Q Abell, Corwin & Olowin (1989); and 4|b Stocke et al. (1991). 


mately 10 and 30 per cent, although some of this variation is 
due to an apparent trend for increasing baryon fraction with 
radius, as shown in Fig. |^. A linear regression to the data in 
this diagram (shown by the dashed line) indicates that the 
baryon fraction may be consistent with Qb.max < 0.06 only 
at the very centre, but the mean value of the data points is 
much higher than the standard primordial nucleosynthesis 
value. We note that Fig. does not account for errors in the 
gravitational potential from the velocity dispersion, but we 
shall address this point in Section 3.2. 


The uncertainty from the core radius, and other param¬ 
eters, on the determination of the baryon fractions has been 
assessed using the Abell 478 data as a control data set. The 
results of these tests, which will be discussed in the following 
section and shown in Table indicate that the gravitational 
potential of the cluster provides the main uncertainty in the 
baryon fraction determinations. 


3 BARYON FRACTION UNCERTAINTIES 

As the deprojection estimates of the cluster baryon fraction 
are given by Mgas /Mgrav , we have estimated the susceptibil¬ 
ity of the deprojection results to uncertainties in Mgas and 
Mgrav resulting from changes in the input parameters for an 
individual cluster. We have also estimated the uncertainties 
in the baryon fraction due to Mgrav using the observational 
errors in the X-ray temperatures. 


3.1 Gas mass uncertainties 

The deprojection method produces gas mass estimates that 
are statistically very well determined. We assume that the 
emission in the outer regions of clusters arises from thermal 
emission rather than non-thermal processes, as there is no 
evidence for significant non-thermal emission at large radii 
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Table 2. Results 


No. 

Cluster 

Ro 

(Mpc) 

dR 

(Mpc) 

Mass (10^"^ 
Gas 

Mo) 

Grav 

Mass Ratio (Mgas /Mgi-av %) 
(R < 1 Mpc) (R < Ro) 

1. 

A85 

1.415 

0.101 

0.87 ±0.06 

4.64 

17.3 ± 1.1 

18.8 ± 1.3 

2. 

A401 

1.265 

0.141 

1.32 ±0.07 

10.1 

12.8 ±0.4 

13.0 ±0.7 

3. 

A478 

1.951 

0.163 

2.38 ±0.21 

9.28 

23.1 ±0.9 

25.6 ± 2.2 

4. 

A545 

1.815 

0.259 

1.91 ±0.25 

10.6 

17.1 ± 1.6 

18.1 ± 2.4 

5. 

A644 

1.198 

0.133 

0.95 ±0.06 

9.06 

10.6 ±0.6 

10.5 ± 0.6 

6. 

A665 

2.376 

0.297 

4.37 ±0.46 

22.1 

18.1 ± 1.0 

19.8 ±2.1 

7. 

A1413 

1.715 

0.245 

1.83 ±0.23 

15.9 

10.8 ± 1.1 

11.5 ± 1.4 

8. 

A1650 

1.090 

0.156 

0.75 ± 0.08 

6.37 

11.8 ± 1.2 

11.8 ± 1.2 

9. 

A1689 

1.481 

0.296 

2.12 ±0.16 

15.5 

13.0 ±0.5 

13.7 ± 1.0 

10. 

A1763 

1.823 

0.304 

2.61 ±0.22 

13.2 

17.6 ± 1.2 

19.8 ± 1.7 

11. 

A1795 

1.426 

0.119 

1.13 ±0.08 

5.49 

18.7 ± 1.1 

20.6 ± 1.5 

12. 

A2009 

1.297 

0.259 

1.44 ±0.10 

10.6 

13.4 ±0.5 

13.6 ± 0.9 

13. 

A2029 

1.291 

0.143 

1.26 ±0.11 

10.3 

11.9 ±0.8 

12.3 ± 1.1 

14. 

A2142 

1.931 

0.276 

2.84 ±0.15 

20.1 

11.9 ±0.3 

14.1 ± 0.6 

15. 

A2163 

2.264 

0.323 

5.46 ± 0.49 

32.5 

14.4 ± 1.0 

16.8 ± 1.5 

16. 

A2319 

1.402 

0.108 

1.73 ±0.12 

14.2 

11.9 ±0.8 

12.2 ± 0.8 

17. 

A3186 

1.508 

0.188 

1.76 ±0.23 

9.50 

15.6 ± 1.9 

18.5 ± 2.4 

18. 

A3266 

1.420 

0.114 

1.42 ±0.07 

9.07 

15.4 ±0.4 

15.7 ±0.8 

19. 

A3888 

1.118 

0.279 

1.20 ±0.15 

8.66 

13.9 ± 1.8 

13.9 ± 1.8 


This table summarizes the deprojection results, where Rq is the outer radius of the deprojection, dR is the bin size. The gas and 
gravitational results are plotted against Ro in Fig. ^ The baryon fractions within 1 Mpc and the total region of each deprojection are 
given in the last two columns. Note the observational errors in the velocity dispersion are not available for all the clusters, and so are 
not quoted. The uncertainty in Mgas and Mgas /A^grav are Icr standard deviation values, resulting from the statistical uncertainty in the 
X-ray data. 


from the radio waveband. The main uncertainty in the gas 
masses arises from the intrinsic X-ray luminosity of a clus¬ 
ter, i. e. through the estimate of the distance to the cluster, 
intervening absorption, spherical symmetry, and the effect 
of clumping in the intracluster gas. All these points are ad¬ 
dressed below. 

The effect of ellipticity in the cluster X-ray emission has 
been investigated by D. White et al. (1994) in their analysis 
of ROSAT HRI data on A478. They found that the ellip¬ 
ticity of (1 — h/a) = 0.2 in the X-ray emission produced an 
average (and la ) value of Mgas = (4.6 ± 0.5) x 10^^ Mq 
( within 0.5 Mpc) from the deprojection of four sectors, as 
compared to Mgas = (4.8 ± 0.2) x 10^® Mq from an az¬ 
imuthal average. Thus, within errors the effect of the spher¬ 
ical symmetry assumption is negligible. We also note that, 
although a cluster may appear spherically symmetric in pro¬ 
jection, it may be extended in the line of sight. However, 
for a constant luminosity Lx oc Mgas ^/V the volume V 
would have to be increased by a factor of 16 to eliminate 
baryon over-densities of 4. Similarly, the accuracy of the 
background subtraction, which affects the luminosity esti¬ 
mate, would have to be wrong by a factor of 16 to reduce a 
baryon fraction of 25 per cent to the universal value of < 6 
per cent. We therefore do not consider spherical asymme¬ 
tries, either tangential or elongation along the line of sight, 
or background subtraction uncertainties, to be important ef¬ 
fects in the baryon overdensities in clusters, especially if the 


Figure 1. 

This diagram shows that differing mass fraction profiles obtained 
with differing core radii (0.2, 0.5 and 1.0 Mpc) for the gravita¬ 
tional mass distribution. This example is for the Abell 478 data, 
where the core radius used to give a fiat temperature profile 
was 0.2 Mpc. This is also approximately the core radius deter¬ 
mined from a comparison of a deprojection and spectral analysis 
of ROSAT PSPC data (Allen et al. 1993). We note that outside 
the core region of each potential the mass fraction profiles are ap¬ 
proximately fiat, and more importantly, tend to the same result. 

Figure 2. 

This diagram shows the baryon fraction (Mgas /Mgrav ) at the 
outer radius of each deprojection. The dashed line show a best¬ 
fitting linear function of Mgas /Mgrav = 0.0579 -I- 0.0556R. This 
is clearly inconsistent with the standard nucleosynthesis value of 
< 6 per cent, indicating by the dot-dashed line. Note the dashed 
line also shows an increase in the baryon fraction with radius. 
Observational errors on Mgrav are not included in t his plot but 
the effect on Mgas /Mgrav is estimated in Section |3.2| from Fig. 

baryon overdensities are shown to be common in a statistical 
sample of clusters such as ours. 

The uncertainty in the gas masses from the distance 
is obviously dependent on cosmological parameters and the 
cluster redshift (we have adopted a Hubble constant of 
Ho = 50L5okms“^ Mpc”*^ in the general analysis). The 
expected dependences of the masses on Ho are approxi¬ 
mately Mgas oc h^o^^ ^ constant radial density profile 

in the cluster), Mgrav oc h'^Q (at radii outside the core of an 
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Table 3. Test Parameters 


Test 


Cosmology 

JV„ 

t-k 


dM/dR 

cr 

^core 

Po 

■^gas 

■^grav 

■^gas / -^grav 

No. 

Ho 

90 

z 

(10^^ cm“2 ) 

(keV) 

G-C 

(Mq kpc“l) 

(kms“l ) 

( Mpc) 

(10^ Kcm“3 ) 

(lOl^Mg) 

0 

© 

(xioo%) 

0. 

50 

0.0 

0.0881 

1.36 

6.8 

ISO-ISO 

N/A 

904 

0.20 

1.5 

2.38 ± 0.21 

9.28 

25.6 ± 2.2 

1. 

100 








0.10 


0.84 ± 0.08 

4.65 

18.2 ± 1.6 

2. 


0.5 









2.30 ± 0.20 

9.11 

25.2 ± 2.2 

3. 



0.0890 








2,44 ± 0.21 

9.35 

26.1 ± 2.3 

4. 



0.0872 








2.32 ± 0.20 

9.21 

25.2 ± 2.2 

5. 




2.50 







2.59 ± 0.23 

9.28 

28.0 ± 2.5 

6. 





7.9 





2.0 

2.38 ± 0.20 

9.28 

25.7 ± 2.2 

7. 





5.8 





1.0 

2.38 ± 0.21 

9.28 

25.6 ± 2.3 

8. 






KNG-KNG 




2.0 

2.38 ± 0.21 

7.08 

33.6 ± 2.9 

9. 






NO-ISO 





2.38 ± 0.21 

8.10 

29.4 ± 2.6 

10. 






NO-LM 

5.0 

N/A 

N/A 


2.38 ± 0.21 

10.2 

23.3 ± 2.0 

11. 








1165 

0.50 

1.0 

2.39 ± 0.21 

17.1 

13.9 ± 1.2 

12. 








764 

0.15 

3.0 

2.39 ± 0.20 

6.67 

35.8 ± 3.0 


This table summarizes the effects of uncertainties in various input parameters used in the deprojection analysis on the results (shown in 
the last three columns). The tests have been applied to the A478 data, and the variations should be compared with the standard results 
shown in the first row (test number 0). The largest reduction in the mass ratio is produced by lowering the velocity dispersion to the Icr 
lower limit given by Zabludoff, Huchra & Geller (1990). The parameter labeled (p indicates the galaxy-cluster combined potential used; 
ISO indicates a true isothermal potential, KNG a King Law potential, NO a null contribution, and LM indicates a linear mass model. 
The numbers for the gravitational potentials are: a for the velocity dispersion of the cluster and rcore for the core radius, or dM/dR for 
the amount of mass in the linear mass model. Pq is the pressure used at Rq to obtain the correct deprojected temperature profile (in 
conjunction with the core radius where applicable). N/A indicates the entry was not applicable to the potential used in that test. 


isothermal sphere), and therefore the baryon fraction should 
change as Mgas /Mgrav oc ■ However, we have found 

that a deprojection with a different Hubble constant requires 
a gravitational potential with a proportionately smaller core 
radius to obtain the same temperature profile as that ob¬ 
tained with a smaller Hubble constant. Test number 1 of 
Table ^ shows that with hso = 2 the change in Mgrav is 
in agreement with that expected for a cluster that is half 
as distant and with a core-radius half as large. The corre¬ 
sponding change in Mgas is less than the expected value of 
0.42 X 10^^ M©, because the required change in core radius, 
for a flat temperature profile, results in a larger X-ray lumi¬ 
nosity and gas content in the central regions of the cluster. 
Therefore, because the changes in the Hubble constant force 
a recalibration of the deprojection results, the Hubble con¬ 
stant uncertainties lead to a smaller changes in the baryon 
mass fraction than would be expected. The uncertainty in 
go and the redshift of a cluster produce comparatively small 
changes, as shown in test 2 for go = or tests 3 and 4 for 
the statistical uncertainties in the redshift of A478. 

Although the Hubble constant provides the greatest 
uncertainty in the gas mass determinations, it does not 
eliminate the large baryon over-densities. An unreasonably 
small Hubble constant would be required to reduce them to 
the standard primordial nucleosynthesis values, because this 
also depends on Ho as Hb < O-Odh^Q . However, as Steigman 
(1987, 1989) has noted, a more useful limit may be obtained 
by requiring that the gas mass does not exceed the total 
mass of the cluster. Assuming that Mgas/Mgrav oc 
then we obtain a lower limit on the Hubble constant of 
Ho = 22 kms“^ Mpc“^. 

The gas mass determinations also depend on the esti¬ 
mate of the absorption of X-rays emitted from the cluster. 
We have already noted that intrinsic absorption may occur 


Figure 3. 

These plots show: (a) the error in the gas mass estimate, (b) 
the emission-weighted temperature, when the X-ray emission is 
assumed to be from a single-phase medium but there are actu¬ 
ally two phases. The single-phase temperature is assumed to be 
kTref = 7 keV. The main-phase temperature is kTi = 7 keV (with 
an abundance of Zi = 0.4 Zq), and the secondary-phase temper¬ 
ature is varied between kT 2 = (0.01 — 10) X kT^ef. The separate 
lines are for volume filling factors of the secondary phase of V 2 - 
0.0 — solid (flat), 0.01 - dash, 0.05 - dash-dot, 0.10 - dot, 0.70 - 
dash-dot-dot-dot, 0.5 — solid. 

Figure 4. 

These plots are similar to those in Fig. ^ where the single-phase 
temperature is assumed to be kT,.ef = 7keV, but here the main- 
phase temperature is actually kTi = 15 keV and the secondary- 
phase abundance is Z 2 = 2 Zq . 

in some clusters, but we have used Galactic column den¬ 
sities determined from 21 cm measurements throughout to 
give a consistent sample of column density determinations. 
As these represent minimum estimates for the total column 
densities, we have tested for the effect of excess absorption 
on the baryon fraction. We expect that the baryon fraction 
should increase with the intrinsic luminosity, and therefore 
gas masses will be larger rather than smaller. A478 provides 
an ideal example, as the excess absorption in this cluster 
has been well studied (D. White et al. 1991b, Johnstone et 
al. 1992, Allen et al. 1993). In test 5 of Table ^ we show 
that an addition of 1.1 x 10^^ cm“^ above the Stark et al. 
(to the value determined by Allen et al. from their spectral 
fits of the ROSAT PSPC data on A478 of 2.5 x 10^*^ cm“^ ) 
produces approximately a 10 per cent increase in the gas 
mass. 

One further point in the determination of gas masses 
from X-ray data that we discuss is the effect of clumping 
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in the intracluster gas. We have estimated the error in the 
determinations of the gas mass that could arise when the 
gas is assumed to be a single-phase medium, but in actual¬ 
ity the gas is multiphase. Two phases have been considered, 
and the combined emission is forced to produce a fixed total 
number of counts To. 4-4 kev in a waveband from 0.4 — 4 keV 
{i.e. a top-hat approximation to the response of the IPC). 
We then select a reference temperature for the single phase 
estimation and compare this mass with the mass that we 
would estimate if the gas had two-phases with different tem¬ 
peratures and volume filling-factors. The masses in the two 
phases are given by solving the following equation assuming 
pressure equilibrium between the two phases: 

/*4 keV 


-^0.4-4 keV Cx: TliVi 


A(kri 


/0.4keV 

/*4keV 


nlV2 


E 
A(kr2 


dE + 


■dE 


( 2 ) 


where the subscript number refers to the two phases, n is 
the electron number density, kT is the temperature variable, 
V is the volume fraction, and E is photon energy in the 
integral that evaluates the total number of counts from the 
cooling function A in the specified waveband. We note that 
equation ^ takes no account of absorption or the effect of 
cluster redshifts. Note, we also assume pressure equilibrium 
between the two phases, otherwise a mechanism is required 
to prevent the cooler gas from expanding and mixing into 
the hotter gas. 


In Fig. ^(a) we show the error in the gas mass determi¬ 
nation when a single phase of temperature kTref = 7keV is 
assumed. The lines show the mass error when there is one 
component of temperature kTi = 7keV and a secondary 
phase which is varied between kr 2 = (0.01 — 10) x kT^ef. 
The different lines show the mass error for volume fractions 
of the secondary phase, V 2 = 0 — 0.5. It can be seen that 
a baryon over-density of a factor of 2 could be eliminated 
by over-estimates in the gas mass determination, if the sec¬ 
ondary phase filled less than approximately 40 per cent of 
the total volume, and had a temperature between approxi¬ 
mately 0.8 and 1 keV (depending on the exact value of V 2 ). 
However, from Fig. ib), we can see that the corresponding 
emission-weighted temperature from the combined medium 
could only be as high as approximately 1.5 keV, irrespective 
of V 2 , so that it is unlikely such errors in the gas mass could 
be made as the temperature was assumed to be kTi = 7 keV. 
Observational uncertainties would usually rule out such a 
large discrepancy. 


From a slightly different perspective, one can ask if 
we can obtain sufficient gas mass overestimates when the 
emission-weighted temperature from the combined emission 
is close to that expected from a single-phase gas. In Fig. ^ we 
show the results when the gas is thought to have a single¬ 
phase temperature of kTref = 7keV, but there is actually 
a component at kTi = 15keV (of the same abundance 
of Zi = 0.4 Zq) and a second component, again between 
kr2 = (0.01 — 10) X kTref (this time with an abundance of 


Figure 5. 

This diagram shows the different gravitational mass distributions. 
The standard deprojection results employ the true isothermal po¬ 
tentials (solid line). We note that the King law underestimates 
the mass at outside 8 to 10 core-radii (which is =0.2 Mpc in this 
example). 


Z 2 = 2.0 Zq). Very large overestimates can be produced, 
but a factor of two overestimation is not obtained unless 
the emission-weighted temperature is allowed to be as low 
as approximately 5keV (for V 2 = 0.01). In this case the av¬ 
erage abundance would be about 0.5 — 0.6 Z©, which is not 
unreasonable compared to the Z^^f = Zi = 0.4 that would 
be assumed, and the fraction of mass in the cooler phase is 
approximately 10 per cent (the luminosity contribution is 
about be 70 per cent). 

We can apply this scenario of significantly different tem¬ 
perature phases to a cluster of a similar emission weighted 
temperature. A1763 has an emission-weighted temperature 
of kT ^ 7keV, and a deprojected gas mass of 2.6 x 10^'* Mq 
( within 1.8Mpc radius). Therefore, from our example, we 
would expect 1.3 x 10^^ M© in a cooler phase to produce a 50 
per cent overestimate of the gas mass. This amount of cooler 
gas cannot be contained within the interstellar medium of 
giant elliptical galaxies (which have the required tempera¬ 
ture of approximately IkeV), as the mass in the cool gas 
is equivalent to approximately a thousand giant elliptical 
galaxies, which is clearly unreasonable. Thus the majority of 
the cooler gas would have to be in the intracluster medium, 
isolated from the destructive processes of the hotter phase 
by magnetic fields. The problem with this scenario is that 
observations already appear to rule out variations of more 
than a factor of two in the intracluster gas, as we discuss 
below. 

In summary, large gas mass overestimations can occur 
when there is significant amounts of cooler gas at kT^l keV 
with line emission which enables the same count emissiv- 
ity to be produced by a smaller mass of gas. As the effect 
is due to the lines, the abundance of the intracluster gas 
also influences the possibility of mass determination errors. 
However, the emission-weighted temperature also decreases 
with abundance, as most of the emission comes from the 
cooler phase, and the resulting effect of abundance varia¬ 
tions is that the gas mass over-estimates are very nearly 
constant for a given range of the emission-weighted temper¬ 
ature. We note that, in a spectral analysis a contribution 
from a cool phase should be easily discernible, however our 
deprojection analysis is a broad-band analysis and cannot 
discriminate between combined spectra of various tempera¬ 
ture which produce similar count emissivities. 

Although we cannot rule out such disparate tempera¬ 
tures from our imaging analysis, a spectral analysis of Ginga 
and EXOSAT data on the Perseus cluster (Allen et al. 1992) 
only allows temperature variations of a factor of approxi¬ 
mately two. Also, a spectral analysis of the A478 cluster 
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Figure 6. 

These plots show how we have estimated the effect of uncertain¬ 
ties in the gravitational potential using the errors in the observed 
X-ray temperatures (from 13 clusters where the temperature er¬ 
rors are measured, and symmetric to within 2keV). The uncer¬ 
tainty in the gravitational mass has been estimated by propa¬ 
gating the (Icr) errors in the observed X-ray temperature to the 
baryon fraction at 1 Mpc, as shown in (a). Assuming that the er¬ 
rors are symmetric and Gaussian we have then determined the 
cumulative probability, as shown in (b), from which we estimate 
that the cluster baryon fraction at 1 Mpc has a median value of 
13.8 per cent, with 5th and 95th percentile limits of 10.0 and 22.3 
per cent. 

out to 2 Mpc (Allen et al. 1993) indicates that a 1 keV com¬ 
ponent cannot be significant in this cluster, as the best-fit 
emission-weighted temperature is consistent with the broad- 
beam value (6.8 keV), and is above 4keV at the 90 per cent 
confidence level. Only in the central regions of the cool¬ 
ing flow, and between 1 — 2 Mpc is the lower-limit around 
1 keV (but the best fit is around 4keV). Thus, within 1 Mpc 
where the temperature is well constrained and there is still 
a baryon over-density problem, the results indicate that a 
cool component is not significant. We expect ASCA to be 
able to rule out such variations to much larger radii. 

One further point is that the baryon overdensities are 
common to the whole sample and do not appear to be de¬ 
pendent on the Galactic column density. If clumping were 
responsible for gas mass overestimates then we would have 
expected clusters such as A478, which have large Galactic 
column densities, to have smaller than average baryon over¬ 
densities because we would see little of the sub-1 keV emis¬ 
sion would be responsible for the overestimations. 

From our investigations into the required conditions for 
significant gas mass overestimations, and spectral observa¬ 
tions of specific clusters, we conclude that clumping cannot 
explain the baryon overdensities in clusters. 

3.2 Gravitational mass uncertainties 

We have shown that the gas mass uncertainties are unlikely 
to reduce the cluster baryon fractions to the 6 per cent up¬ 
per limit obtained from standard primordial nucleosynthe¬ 
sis. However, the gravitational potential is the most uncer¬ 
tain component in the calculation and we now discuss its 
uncertainties. The deprojection results are changed by al¬ 
tering the gravitational potential to give a temperature that 
is consistent with the observed broad-beam values. In tests 
6 and 7 of Table we can see that the statistical uncertain¬ 
ties in the temperature (for A478) produce comparatively 
small changes in the results, so that individual baryon frac¬ 
tion uncertainties will be probably dominated by the form 
of the potential that is chosen to obtain this temperature, 
rather than errors in this temperature determination. 

We have investigated the effect of changes in the grav¬ 
itational potential, e.g. for a King-law density distribu¬ 


tion, the gravitational contribution from the central galaxy, 
changing the form of the gravitational potential, and the sta¬ 
tistical uncertainties in the optical velocity dispersion of the 
cluster. In the first case, test 8 shows that replacing both the 
galaxy and cluster potentials with King-law distributions in¬ 
creases the baryon fraction estimate. The reason for this is 
shown in Fig. ^ where we have plotted several different mass 
distributions (appropriate for Abell 478, i.e. a velocity dis¬ 
persion of 904 km s”*^ and the core-radius which we used of 
0.2 Mpc). The King approximation provides a good descrip¬ 
tion for the mass distribution within 10 core-radii (< 2 Mpc), 
but outside this region the King law clearly underestimates 
that total gravitational mass compared to the true isother¬ 
mal potential. Although we have no particular reason to 
believe the cluster should follow a true isothermal potential 
at large radii (especially if the cluster in not relaxed), we use 
the true isothermal potential to provide conservative baryon 
fraction estimates. In test 9 we show that when the mass of 
the central galaxy is neglected, a larger baryon fraction is 
estimated for the cluster. When the remaining cluster po¬ 
tential is changed to a linear mass distribution (test 10), 
similar results to the standard result are obtained. 

The major change in the gravitational mass estimates, 
and therefore the baryon fraction, actually arises from the 
uncertainties in the velocity dispersion, as shown in tests 
11 and 12. The statistical uncertainties (for A478; Zablud- 
off, Huchra & Geller 1990) indicate that the cluster baryon 
fraction can be reduced from 26 per cent to 14 per cent 
when the (In) upper limit on the velocity dispersion is used 
(-1-261 kms“^), but it is then very difficult to obtain a fiat 
temperature profile, and the baryon fraction has still not 
been reduced to less than 6 per cent. To reduce all the baryon 
overdensities to < 6 per cent would require that we use high 
velocity-dispersions for all the clusters, and would then pro¬ 
duce unsatisfactory temperature profiles. This seems to be 
an unlikely solution to the baryon overdensity problem, espe¬ 
cially as optical velocity dispersions are, if anything, usually 
overestimated due to substructure in clusters. 

We have attempted to estimate the uncertainty in the 
results due to the gravitational potential mass using the er¬ 
ror in the reference (observed) temperatures. To estimate 
the effect of uncertainties in the velocity dispersion we have 
plotted the baryon fraction within a consistent radius of 
1 Mpc (see Table against the observational X-ray tem¬ 
perature (see Table ^). We have only included those data 
where the X-ray temperature has been measured and its 
uncertainties (the standard deviation errors) are reasonably 
symmetric (i.e. the positive and negative errors are similar 
within 2keV). This eliminates A545, A1689, A1763, A2009, 
A3186 and A3888. Using this refined sample of 13 clusters 
we have propagated the uncertainty in the temperature onto 
the uncertainty in the baryon fraction, as shown in Fig. ^(a). 
(Note, Mgas is not significantly affected by uncertainties in 
the X-ray temperature.) To then estimate confidence limits 
on the baryon fraction with 1 Mpc we have treated the errors 
as Gaussian, and determined the cumulative probability as 
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Figure 7. 

The gas masses (squares) at the outer radius of the deprojection 
of each cluster are plotted together with the total gravitational 
masses (triangles). The solid line shows a fits to the gas masses, 
Mgas = 6.7X 10^^-^^Mpc’ which predicts gravitational masses [dot- 
dash line) of Mgrav (S^b — 0.06) = 1.1 x 10^^if flb/Ho = 
0.06. The actual gravitational masses in the deprojection results 
are fit with Mgi-av ' = 4.7 x [dashed line) if the same 

index as the gas mass is used, or Afgrav " = 5-3 x if the 

power law has a free-fit index [dotted line). (Errors on the data 
points are 1 standard deviation.) 

a function of baryon fraction, as shown in Fig. |(b) . The di¬ 
agram clearly shows that although there is a wide variation 
in the baryon fraction, it is very unlikely {i.e. a probability 
of 10“"^) that at least one cluster from a similar sample has 
a baryon fraction at 1 Mpc of fib,max < 0.06. We estimate 
that the median baryon fraction at 1 Mpc is 13.8 per cent 
with 5th and 95th percentile confidence limits of 10.0 and 
22.3 per cent. 


4 DISCUSSION 

The results, as shown in Fig. from our deprojection analy¬ 
sis of 19 clusters of galaxies indicate that the baryon fraction 
in clusters is inconsistent with the mean baryon fraction of 
the Universe predicted from standard primordial nucleosyn¬ 
thesis calculations, if flo = 1, as first noted by S. White & 
Frenk (1991) for the Coma cluster. The diagram also shows 
a trend for increasing baryon fractions with radius, and in¬ 
dicates that the cluster baryon fraction could be consistent 
with the universal value of < 6 per cent at the very centre 
of clusters, but not further out. 

Our determinations of the baryon content in clusters are 
not compromised by the uncertainties in our analysis. The 
gas masses are extremely well determined, and overestimates 
due to clumping appear unlikely to be able to simultane¬ 
ously reduce the baryon fractions significantly and produce 
observationally consistent emission-weighted temperatures. 
The uncertainty in Hq would require an unreasonably small 
Hq to reduce the cluster baryon fractions to < 6 per cent 
(flo = 1), and a lower limit of Hq = 22 kms“^ Mpc“^ is ob¬ 
tained by allowing all the mass of the clusters to be in gas 
(see also Steigman 1987, 1989). Possible excess absorption 
in clusters only increases the gas mass estimates. The main 
uncertainty in the cluster baryon fractions probably arises 
from the uncertainties in the total gravitational mass, which 
are dominated by the cluster velocity dispersion values. The 
optical determinations of all the velocity dispersions would 
be have to be underestimated, which is somewhat contrived, 
and is also contrary to overestimates expected from optical 
determinations if the clusters have undetected substructure. 

Since our results indicate that baryon fractions at 1 Mpc 
are typically 10 — 20 per cent in clusters, then the simplest 
solution to the conflict with standard primordial nucleosyn¬ 


thesis may be that f2o^0.3. As there is evidence for flo = 1 
(see the Introduction) we shall first discuss the implications 
that arise from assuming standard primordial nucleosynthe¬ 
sis when flo = 1. (Note we ignore the fact that go = ^ when 
flo = 1, whereas our results are for go = 0. Test 2 in Table ^ 
indicates that go has little affect on the results.) 

Using the results given in Table ^ we have plotted, in 
Fig. 0. the gas and gravitational masses against the max¬ 
imum radius of each deprojection. We note that if all the 
cluster deprojections are extended to the surface brightness 
of the background, then we would expect the gas masses at 
the maximum radii to follow an Rq dependence, and indeed 
fitting a power-law function to the gas masses at the maxi¬ 
mum radius indicates that the index is 2.2 with 90 per cents 
confidence limits of ±0.17. We have therefore obtained gas 
masses approaching the maximum detectable radii for these 
data. Forcing an index of 2, and fitting a power law to the 
gas masses gives the fit of Mgas = 6.7 x lO^^i^^pc shown 
by the solid line. The corresponding total masses used in 
the deprojection analysis, fitted with the same radial de¬ 
pendence, gives Mgrav ' = 4.7 X 10^‘^i?Mpc, shown as the 
dashed line. If the gravitational masses are fitted with the 
radial dependence as a free parameter, then we find that 
Mgrav " = 5.3 X 10^'*7?Mpci showu as the dotted line. This 
again indicates that the baryon fraction increases with ra¬ 
dius, as already found in Fig. If we then assume that 
flb/flo = 0.06 in clusters, and return to the same radial de¬ 
pendence as the gas masses, then the expected total gravita¬ 
tional mass is given by Mgrav (fib < 0.06) = 1.1 x 10^®J?Mpci 
shown as the dot-dash line. 

From this we can see that, if we use the mean baryonic 
fraction of < 6 per cent to predict the total gravitational 
masses from gas masses, then we overpredict masses with 
respect to the virial values, e.g. for A665 the predicted mass 
is approximately 5.7x 10^® Mq within 2.4 Mpc. This is larger 
than considered from current theories of cluster formation, 
which would give a total mass of 2.8 x lO*^® M© for A665 
[from equation kr/(4keV) = (M/10^® M©)^/® in Henry et 
al. 1992]. This total mass is more in line with that expected 
for a very much hotter cluster, such as A2163 at 13.9 keV. 
We also note that current theories of the formation of large- 
scale structure and cluster of galaxies may have problems ex¬ 
plaining the apparently common occurrence of large baryon 
overdensities. For the median and 5th and 95th percentile 
confidence limits that we have placed on the baryon fraction 
within 1 Mpc, the overdensity is probably at least 2flb. 

If clusters are truly overdense in baryons then, as high¬ 
lighted by Fabian (1991) from the Shapley Supercluster 
data, then how are the extra baryons accumulated from 
the surrounding volume at the maximum mean density of 
6 per cent for the Universe? In A665 the gas mass is ap¬ 
proximately 5.4 X 10^"* M© within 2.4 Mpc, and therefore 
the size of the region with the equivalent mass of baryons 
for a Universe of density fib < 0.6 is 31 Mpc — a factor of 
13 in radius, or greater than 2 x 10® in volume. The require- 
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ment of such large regions, to provide a source of baryons 
for such overdensities, may be too large to enable the con¬ 
centration to occnr within a Hubble time and would rule out 
self-gravitational accumulation of baryons as a valid mech¬ 
anism to concentrate the baryons. This is the real problem 
of baryon overdensities in clusters, as it is independent of 
the uncertainties in the gravitational mass estimates in this 
analysis. However, even if we assume that sufficient baryons 
can be accumulated within the cluster, we still have to ex¬ 
plain how the baryons appear be concentrated at the centre 
of a cluster with respect to the overall dark matter distribu¬ 
tion. In Fig. Ijwe show a schematic diagram of the gas and 
gravitational mass distributions that could give rise to large 
baryon fractions within the central ~ 3 Mpc, decreasing to 
a baryon fraction consistent with the universal average at 
larger radius. [We note, with the mass fractions determined 
from the deprojection results and the ft values for clusters 
(Forman & Jones 1984), both indicate an increase of the gas 
to gravitational mass fraction increases with radius, over the 
observed regions of clusters, i.e. ^3 Mpc]. 

We can envisage two ways to create the distribution 
shown in Fig. ^ - through evolution or an uneven distri¬ 
bution of baryonic and non-baryonic material in the early 
Universe {z ^ 5). First, an evolutionary process, which 
may produce a central concentration of gas surrounded by 
an ‘extended halo’ of dark matter, from the infall process 
which forms the clusters and/or the subsequent infall of sub¬ 
clusters. For example, if a gas-rich cluster fell into a larger 
cluster the gaseous component would be stripped from it in 
the dense central regions of the larger cluster, in a manner 
similar to the ram-pressure stripping of the hot gas from 
the elliptical galaxy M86 in the Virgo cluster {e.g. D. White 
et al. 1991a), while the collisionless dark matter would pass 
through unhindered to the other side of the cluster. This 
process would produce an atmosphere of gas which would 
be slightly more extended than the virial core of the clus¬ 
ter, due to shock heating, surrounded at larger radius by a 
halo of dark matter. This dark matter, if bound, may re¬ 
main at large radius for a relatively large period of time 
before falling again towards the core of the cluster. Thus, 
within the framework of hierarchical merging, infalling sub¬ 
clusters may produce significant amount of dark matter at 
large cluster radii. This scenario requires that clusters are 
more massive, approaching 10^® M©, than generally consid¬ 
ered in current theories of cluster formation, and would re¬ 
sult in large peculiar velocities around massive clusters of 
galaxies. Other methods in which the dark matter could be 
distributed on larger scales rely on different clustering prop¬ 
erties of the dark matter, e.g. if the Universe is composed 
of a mixture of mostly hot with some cold dark matter, or 
if A is non-zero. 

Alternatively, if the central concentration of baryons 
with respect to the dark matter does not occur in the evolu¬ 
tionary scenario, then the gas needs to be distributed differ¬ 
ently before the formation of clusters. However, as the gas 
is more concentrated than the gravitational matter, gravita¬ 


tional effects cannot have been responsible, and the baryons 
must have been pushed together to form regions of higher 
density. This could have happened if there was a population 
of active quasars with strong winds or radiation pressure 
which produced voids in the early Universe before cluster 
formation. The baryonic material would have been forced 
together at the interface between voids, at the sites of clus¬ 
ter formations, while the dark matter would have been less 
compressed. Clusters would then have inherited the distribu¬ 
tions of baryonic and non-baryonic material. This scenario 
leads to the prediction that there should be a population of 
objects at the centre of voids. 

None of the above solutions for the baryon over¬ 
densities resulting from standard primordial nucleosynthe¬ 
sis are very elegant or without problems. Perhaps the most 
damning fact is that it appears extremely difficult to accu¬ 
mulate enough baryons from a region with a baryon density 
of at most 6 per cent to provide the overdensity seen to be 
common in our sample. As there is still evidence for Hq = 1 
on large scales, e.g. from studies of the structure in clusters 
(Richstone, Loeb & Turner 1992), and the POTENT analy¬ 
sis of /7 ?j 45galaxies (Nusser & Dekel 1993, Dekel et al. 1993, 
Dekel & Rees 1994). We do not appeal to low values of Ho, 
but assume that the dark matter in clusters is spread over 
a larger radius than the baryons. This means that clusters 
are several times more massive than is canonically assumed. 


5 CONCLUSION 

Our deprojection analysis of 19 moderately luminous and 
distant clusters, observed with the Einstein Observatory 
IPC, shows that cluster baryon fractions are all inconsistent 
with the mean value for the Universe of Hb = 0.05±0.01h/)/, 
as calculated according to standard, homogeneous, pri¬ 
mordial nucleosynthesis (Olive et al. 1990, Walker et al. 
1991). The deprojection method produces well-determined 
gas masses, such that the main uncertainty in the gas mass 
lies in the value of the Hubble constant, while the overall 
main uncertainty in the baryon fraction determinations lies 
is in the gravitational masses. However, this also cannot pro¬ 
duce a significant enough effect to reconcile the cluster de¬ 
terminations with the mean value predicted from standard 
primordial nucleosynthesis. We find, at the 5th and 95th per 
cent confidence levels, that the baryon fractions of the clus¬ 
ters, in our refined sample of 13, lie between 10 and 22 per 
cent. ASCA should reduce uncertainties in the gas and grav¬ 
itational mass determinations, by enabling accurate temper¬ 
ature measurements (with adequate spatial resolution) to be 
made, from which the total masses and baryon fractions of 
clusters will be accurately determined. 

As there is still strong evidence that Ho = 1 on large 
scales, we have considered the implications that result from 
conflict between the baryon fractions in clusters and the 
mean baryon fraction prediction from standard primordial 
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Figure 8. 

This schematic figure shows how the observational results, which indicate baryon fractions approaching 30 per cent, may be reconciled 
with the mean baryon fraction for the Universe of < 6 per cent (Ho = 1 and h^Q = 1) for the cluster as a whole. The solid lines are the 
cluster gas and gravitational mass distributions, and the dotted line shows the mass expected within the same volume with a critical 
density of material and the baryon fraction of 6 per cent (flo = !)• The reasoning for the more extended nature of the dark matter is 
given in the main text. 


nucleosynthesis. These solutions, which imply clusters are 
much more massive than generally thought, require halos of 
dark matter outside the main X-ray extent of the cluster. 
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